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HOC  Spectral  Analysis  of  an  Almost  Periodic  Random  Sequence  in 
Noise 


1 .  Introduction 

Consider  a  stationary  random  sequence  (Zt),  t  ■  0,11,12,..., 
given  by  the  equation 


(A. 


COS  W^t 


+  B. 


sin  w^t) 


+  K, 


j=l 

where  the  amplitudes  A^.Bj,  j  *  l,...,p  are  random  variables  and 
Kt  Is  a  random  colored  noise  Independent  of  the  A^,Bj.  The 
problem  addressed  in  this  paper  is  to  determine  p  and 
from  expected  zero-crossing  counts  regardless  of  the  magnitude  of 
the  noise  term  1 t .  This  can  be  done  under  some  fairly  general 
conditions  by  considering  sequences  of  expected  zero-crossing 
counts  obtained  by  repeated  filtering  of  (Zt). 

Because  zero-crossings  portray  the  oscillation  in  (Zt>  the 
problem  is  tantamount  to  determining  p  and  the  «'s  of  extreme¬ 
ly  weak  signals  burled  In  noise  from  the  oscillation  of  the  process 
and  Its  filtered  versions.  The  general  message  of  the  present  work 
is  that  this  is  indeed  possible. 

More  precisely,  we  shall  be  dealing  with  the  so  called  higher 
order  crossings  (HOC) .  Higher  order  crossings  are  zero-crossing 
counts  observed  in  a  process  and  In  its  linearly  filtered  versions. 
The  filtering  operation  may  consist  of  a  single  operation  applied 
once  or  applied  repeatedly  or  It  may  consist  of  a  succession  of 
different  filters.  In  this  regard,  repeated  differencing  and  re¬ 
peated  summation  play  an  Important  role.  The  shifts  or  changes  In 


the  spectral  distribution  resulting  from  the  filtering  operation 
are  captured  very  economically  by  the  higher  order  crossings,  a 
fact  that  led  to  a  methodology  useful  In  discrete  spectrum  analysis 
as  reviewed  in  Kedem  (1986).  The  present  paper  is  a  refinement  of 
the  work  reported  there.  Specifically  we  Investigate  the  conver¬ 
gence  of  sequences  of  expected  HOC,  as  well  as  functions  of  expect¬ 
ed  HOC,  to  the  (i) 1  s  under  various  conditions  on  the  spectrum  of 
{Kt)  ■ 

When  (Z^)  is  Gaussian  with  mean  zero,  the  oscillation 
depicted  by  the  expected  HOC  Is  equivalent  to  knowing  the  correla¬ 
tion  function.  However,  HOC  have  a  somewhat  more  direct  interpre¬ 
tation  In  terms  of  the  spectrum  as  can  be  seen  from  their  spectral 
representation  that  will  be  discussed  below. 

In  the  present  paper,  both  (Zt)  and  (<t)  are  assumed  to 
be  stationary  and  Gaussian.  The  first  part  of  the  paper  stresses 
aspects  of  filter  design  and  direct  convergence  of  HOC  to  discrete 
frequencies,  while  in  the  second  part,  consisting  of  section  4, 
the  convergence  is  achieved  rather  indirectly  by  certain  functions 
of  HOC.  The  main  results  of  this  paper  are  Theorems  3.6  and  4.1 
stated  in  sections  3.3  and  4,  respectively.  In  Theorem  3 . 5  we 
assume  Is  white  noise,  p  ■  1,  and  construct  a  sequence  form 

expected  HOC  that  converges  monotonlcally  to  regardless  of 

the  signal  to  noise  ratio  (that  is,  the  ratio  of  the  standard  devi¬ 
ation  of  the  almost  periodic  component  to  the  standard  deviation 
of  the  noise) .  The  discussion  leading  to  this  result  points  to 
the  difficulties  that  arise  once  noise  is  added  to  the  almost 
periodic  harmonic  component.  The  basic  idea  underlying  section  4 


is  the  motion  of  AQ-intervals .  These  intervals  are  instrumental 

in  determining  p,  «1 . «p  regardless  of  the  type  of  noise  and 

its  magnitude,  provided  its  spectral  density  is  square  integrable. 


l.l.  The  model  and  definition  of  HOC. 

Let  {X^}  represent  the  almost  periodic  harmonic  component 
with  p  terms,  p  < 


X 


t 


eo.  «jt  *  !j 


sin  «jt) . 


Withou-  Aoss  of  generality  assume 


0<t>  <  «  <...<«•>  <«. 

1  2  p 

The  (Aj},(Bj)  are  taken  as  uncorrelated  normal  random  variables 
such  that 

“i  ■  EBj  *  -•  EAi*j  -  EBiBj  ■ 8  if] 

EAaBj  ■  0  for  all  i.j. 

We  assume  that  t)  is  a  stationary  zero  mean  Gausslna  process 

with  an  absolutely  continuous  spectrum  and  spectral  density 

f  {  ,.j  )  ,  -n  <  s  n.  For  each  t,  ?  t  has  a  normal  distribution  with 

o 

mean  0  and  variance  .  It  follows  that 

zt  -  xt  ♦ 

is  a  stationary  Gaussian  process  with  a  mixed  spectrum  whose  spec¬ 
tral  distribution  function  can  be  espressed  as  a  sum 


F(w)  *  Fx(y)  +  F  (w),  -«  <  w  s  n. 


where  Fv  is  a  right  continuous  step  function  with  jumps  of  size 


Then  repeated 


>/«  ■ 


Let  *  be  the  shift  operator  *Z. 


t  t-1 

differencing  can  be  defined  by  the  operator 


n 


k*0 


Similarly,  repeated  summation  le  defined  by 

n 


7"zt  .  <!♦«.)%  -  X(£)Zt-fc- 


k=0 


Let  (Yt),  t»0,*l,...,  be  any  stochastic  process  and  let  r [ •  ] 


be  the  indicator  function.  Then  the  number  of  zero-crossings  in 


Yl . Yn  is  given  by 


N 


N 


D  «  2^r[YA?0]  -  2^t[Y1iO,Y1  ^0]  -  {[Y^O]  -  *[YH*0]. 


1  =  1 


1  =  1 


n,  ,N 


The  number  of  zero  crossings  in  <V  Zt)tml  is  denoted  by  Dn+1 


n„  ,N 


and  the  number  of  zero  crossings  in  (A  is  denoted  by 


. , D .  Then  D  , ,  and  ^,D  are  examples  of  HOC.  In  general, 
n+l  n+i  n+i 


when  the  linear  operation  is  a  filter  with  transfer  function  H, 


the  notation  Du  is  used  to  signify  the  HOC  corresponding  to  H. 
H 


The 


problem  is  to  determine  p  and  Wj, . . . from  expected 


HOC  such  as  (ED  },  (ED)  and  (ED„) . 

n  n  n 


2.  Some  moment  relations 


The  second  order  moments  of  Vnz^,  AnZ^  can  be  expressed 


quite  compactly  by  introducing  the  following  sequences.  Define 

n  n 


V 


"•<>>  -  Z[")  •  bo,n'e>  “  Z(")(A)eos(9 

J-0  J«0 


Then  p(0)  =  1,  p( 1)  =  V j  =  ^0  and  so  (2.1)  and  (2.2)  Imply  the 

existence  of  a  function  p  (x, ,xn,...,x  ;9)  such  that 

m  l  i  m 


(2.3) 


p(m)  =  plV  .©_;*) 

m  l  i  m 


pm(  .D  ,  ,  .  .  .  ,  V  ,0)  ,  m  =»  1 

mi  z  m 


The  proof  of  this  fact  follows  easily  by  solving  (2.1),  (2.2) 
recursively  for  p(k)  starting  from  p(0)  =  1,  and  noting  that 

bn+l(n'w)  =  <_1>n+1'  bn+l(n'0)  *  lm 

Thus 

n 

p ( n+1 )  =  (-l)n+1  ^  <Dn+lak(n'w)  "  bk(n,n)  )P]c(1>1  ,  .  .  •  ,»k;«) 

k“0 

n 

*  X  (n+l1)ak(n'0)  “  bit<n'°)  •  •  • 'ic®10)  ' 

k=0 


and  therefore  for  each  n,  the  sequences  (p ( 1 ) ,p ( 2 ) , . . . ,p (n) } , 

{ ED.  ,  ED.  ,  .  .  .  ,  ED  ),  and  (E  .D,E  _D,  .  .  .  ,E  ,D)  are  equivalent. 

Lemma  2.1.  For  a  zero  mean  stationary  Gaussian  process,  the 
sequences  (EDj),(EjD)  and  (p( n) )  are  equivalent. 

this  equivalence  relation  shows  the  relevance  of  HOC  in  spec¬ 


tral  analysis.  It  should  be  noted  that  the  correlation  function 
may  be  obtained  from  many  different  HOC  sequences,  not  just  those 
obtained  by  differencing  and  summation. 


3 .  Filter  design  and  convergence  of  sequences  of  expected  HOC 


He  examine  the  effect  of  several  different  filters  on  zero¬ 
crossing  counts  when  the  process  (Z^)  consists  of  the  harmonic 
signal  only  and  also  when  It  consists  of  signal  plus  white  noise. 

That  is,  $  *  c.  where  {£.}  Is  white  noise.  In  this  case 

a2 

f (w)  *  2~ ,  -*  <  a  *  n . 

3.1.  A  complex  filter. 

Define  a  process  {Yt}  by 

Yt  ■  (l+ei0S)nZt. 

Note  that  Y  depends  on  n  and  0.  The  transfer  function  and 
squared  gain  are  given,  respectively,  by 

H(X )  *  (l+e1(*"*,)n 
I  H ( 1 )  |  2  «  4ncos2n(^). 


Observe  that  {Y^}  is  complex. 


Yt  ”  ut  +  ivt 


where 


n 

ut =  Z[3cos(kei2t-i 


k=0 
n 


- 1(2) 


sin(k0)Zt_k. 


k*0 


It  follows  that  EIY^.12  and  Re  EY^Yt+1  can  be  written  in  two 
equivalent  forms  as  follows.  First, 


i 


$ 


E|Ytl 


-it 


4  V.2 

2- 

j=l 


n 

I  H(A  )  |2dFx(X)  +  |H{X)|2dF£(A) 

J  —  7T 

»2np)— 2np]]^ 


_n 


COS 


2n(6-X' 

l  2 


-n 


and 


Re(EYtYt+i) 


Re^ 


.n 


e  1X|H(\)|2dFx(\)  +  e  U  | H( X ) | 2dFe ( X ) 


-n 

P 


-n 


4“  V  2 
2~ 


C08 


2n 


COS (U j  ) 


*2  rw 

+  _£  4n| 
2n 


-it 


2nfX-0l 

cos  — g—  cos(X)d\. 


Second , 


n  n 


E|Y, 


Eu?  *  Svl  ■  Z  [")  (k-J)e)R(lc-J) 

j-0  k«0 


n 


^ak(n,^)R(k) 


k=«0 


and 


n  n 


Re(EYtYt+1) 


8utVi  *  Evtvt+i  *  Z  Z(")[E]c°'nk-3,‘,,R'k'3*1) 

J»0  k=0 

n+l 

]^bk(n,©)R{k)  . 


k=0 


n+1 

E  I  Y  I  ^  ^ 

1  t1  £  a.  <n,e)p(k) 

k=0 


£  b.  (n,0)p  (Z>  ,...,»  ;w) 

k=0  _  _ K 

n 

E  a.  (n,0  )p.  (D.  ,  •  .  • 
k=0  K  K 


£  b.  (n,0)p  (  D  ,  .  .  .  ,.  i>;0) 

k=0  _ _ _ 

n 

^  a.  (n,e)p.  (  ©,  .  .  . 0;0) 
k=0  K  K  K 


2 

a  n 


£  2r  2nf0+c<>j').  2nf©-<»)jll  .  .  ,  c  r  2nfX-0') 

i5i  j Lcos  r~2"j +cos  [-rLJJcos(“3>  +  —  J.,cos  [— Jc 


2  2  f  2nf0-HJ.'l  .  2nfd-«i'11  a  v  r.  ^ 

£°ir°  brTcos  h-LJJ+4J  -s2n(^]« 


coa ( x )  dx 


(3.1) 

When  {Zt>  Is  a  purely  harmonic  process,  (3.1)  provides  a  way  for 
determining  the  w^'s.  BY  choosing  an  arbitrary  9  e  [0,n],  (3.1 

will  converge  to  cos(«r)  for  «r  closest  to  9  as  n — ►  «>. 

More  precisely  we  have 

2 

Theorem  3.1.  Assume  <7£  *  0,  9  e  [0,n]  and  suppose 


\9~io  |  <  min  (|0-o>.|,  2n-u  . -9 ) 

r  >r  J  J 


Then 


Re  EYtYt+1 


COS(t>  )  . 

r 


Proof.  Since  cos(x/2)  is  monotone  decreasing  in  [0,n],  the 


t-t,  ..i  <i^  i 


»< 

'5 


({] 

*%| 


II 


IK 


Jl 


condition  of  the  theorem  implies  that 


'(V*]  >  'c«fP)'- 


j  *  r, 


The  assertion  now  follows  from  the  last  expression  in  (3.1). 


Since  cos(x)  is  monotone  in  [0,w],  obtaining  coa(wr)  is 


equivalent  to  obtaining  *>r.  r  -  1 . p.  Important  special  cases 


occur  for  0  «  0 ,n 


Corollary  3.1.  Assume  o  =0.  Then 


<y  ,  n — >  ® 

P 


n  E  D 

n 


w  1  ,  n — ►  ® . 


Proof.  From  (2.1),  (2.2)  and  ( 3 . 1 ) ,  as  n- 


£  bk(n,  w)pk(D1 .... 

V  =  jfrj9 - — — —  — ►  cos  («  ) 

n+1  n  P 


5]  aif^ n* n )^it(® i '  ■  ’  *  '^ic'n ) 
k=0 


2  (n,  °)p^(  j®*  •  •  •  *  j^® ® ) 


cos ) , 


5]  a.  (n,0)p.  (.©,••./ «® ; 0 ) 

k-0  *  K 


and  note  that  cos(x)  is  monotone  for  x  e  [ 0 ,n ]. 

Unfortunately,  the  method  just  outlined  breaks  down  in  the 
presence  of  noise.  To  realize  the  effect  of  noise  define  first 


2n  fe-M 


fn(x:e) 


{¥]• 


•A  :.\<V  V.TvCVt' 


Then  we  have 


Lemma  3.1.  The  sequence  of  probability  densities  {f  (\;0)} 
satisfies 


n 

fn(X  ;0)cos(X  )d\ — >008(0),  0  s  9  s  n. 
-n 


Proof .  Suppose  0  <  9  <  n.  Then  for  every  e  >  0  there  exists  a 
6  >  0  such  that 


a  ■  inf  |  cos  f^r^l  I  >  b  a  sup  |cosf^-^-l|. 

|  X  -0 1  <6  lx-0  I  2c  l2J 


Therefore 


2n 


f  (X  ,9  )dX  s  -2Trbs 
n  _  2n 


0,  n — ►  <». 


X-0(S£ 


26  a 


Similarly  for  9  *  n 


c  s 


=  inf  (cosf^-^ll  ^d*  sup  |cosf^-^l 
n  -6  <  |  X  |  «  n  |  X  |  <n-c  l  2  J 


so  that  again 


,2n 


4  7TT5 


0,  n — mb, 


|  X  |  <n~c 


26  c 


Since  cos(X)  is  continuous  and  since  fn(X;0)  is  for  each  n  a 
probability  density  function  with  parameter  9  it  follows  that 


n 

r 


lim 
n-*»  J 


(X 


;9 ) cos ( X )dX 


cos (9 ) ,  0  <  9  5  n . 


□ 


Theorem  3.2. 


Assume 


>  0. 


Then  for  every  9  €  [o,nJ 


5v 


3 


8 «? 

$ 


a 


I 


j 


’va 


4 


v  b(n,0  ),<>.(!>  ,  .  .  .  ,P.  ;W) 

k=0  K  *  1  K 


cos(0 ) ,  n — m» 


I  ai.(n,0  Jp,*!). . *>,.;«) 

k«0  K 


and  the  same  holds  If  the  t> j  are  replaced  by  the  ^0,  and  n 


by  0 . 


Proof .  When  ■  9  tor  some  j,  the  claim  follows  by  bounded 


convergence  from  (3.1).  Otherwise  define 


b'  £  max  |  cos 
U  jsp 


P) 


and  note  that  the  method  used  In  proving  Lemma  3.1  yields 


(b- i2n 


COS2"[^i]dl 


0,  n — > qd . 


Apply  now  Lemma  3.1  to  the  ratio  of  the  two  Integrals  In  (3.1) 


Corollary  3.2.  When  crf  >  0, 


n+l0  *  1 ' 


or  equivalently 


"E  n+lD  „  ,Bn«  .  .  _ 

in - *  °'  TFT - *  "• 


Proof.  The  proof  follwos  from  (2.1),  (2.2) 


Corollary  3.3.  a*  *  0  if  and  only  If  for  some  9  e  [0,«],  9  e 


(Wj  ....  .<yp) 


J]  b.  (n,^  ) P j, ( P .  ,  •  •  •  » P w ) 
k*0 


cos (0 ) , 


V  a  (n,0)p.  (0  , . . . ,0.  ;») 

k*0 


and  the  same  holds  If  the  V  are  replaced  by  the  jD ,  and  n 
by  0 . 

3.2.  The  Slutsky  filter. 

For  k,m,n  =  0,1,2,...,  put 


[(l-8)m(l+8)nr  Zt,  t  -  0 , ±  1 _ 


Then  { Y^  }  Is  a  real  valued  stationary  Gaussian  process  with 
mean  0.  The  transfer  function  of  this  linear  operation  Is  given 


Hk(X)  *=  (l-e"1X  )km{l+e_iX  )kn 


with  squared  gain 


#i«i2  _k(m+n),,  ..km,.  , .  kn 

|H.  (X)|  =  2  '( 1-cos  X)  (l+cos  X) 

The  squared  gain  Is  symmetric  and  unlaodal  in  [0,w]  with  a  peak 


occurlng  at 


X  »  cos 
c 


•1  fn-nTl 
[n+mj 


Let  D 


( 1c )  N 

H  denote  the  number  of  zero-crossings  in  {Y^  't«i' 


where  m,n  are  fixed  and  are  chosen  In  accordance  with  a  prespeci¬ 
fied  X  . 
c 

Results  similar  to  those  obtained  In  the  previous  discussion 

can  be  obtained  in  terms  of  the  HOC  {D„  },  k  ■  1,2,...,  the  dif- 

Hk 

ference  being  the  fact  that  unlike  the  case  treated  earlier  the 


convergence  to  the  can  be  expressed  directly  In  terms  of 

expected  HOC.  This  Is  readily  seen  from  the  equality 


cos 


1 


« 

ft 


m 

4 

ft  3 


k+1 

N-l 


P  e"U!Hk(X)|2dP(X) 
-« 

7*  2 

J  |Hk(X)  rdF(\) 


o2  n 


Z  a j[ ( 1-coso  .  )m( l+coso  .  )n]  COBW  [  [ ( 1-cosX )m( 1+cosX )nlkcosXdX 

j  =  lJ  _  J  J  J  n  J  o 


£  °2A  (  l-coso  .)m(  1+coso  )n]k+-l  f  [  (1-cosX  )■{  1+cosX  )n]kdX 
j»l  J  J  j  Jo 


(3.2) 


Theorem  3.3.  Assume  o*  ■  0  and  suppose  m,n  are  such  that 

(l-cos  u  )m(l+cos  «  )n  >  max  {{l-cos  u.)m(l+cos  u.)n). 

j*r  3  3 

Then 


Proof.  From  (3.2) 


<v  *-“**• 


cos(«r) 


and  cos(x)  Is  monotone  in  [0,n]. 

2 

Corollary  4.1.  Assume  «  0  and  suppose  m,n  are  such  that 
xc  Is  sufficiently  close  to  ur.  Then 


«  ,  k - M». 

r 


Proof.  In  [0,n]  (l-cos  X)m(l+cos  X)n  is  unimodal  with  the  peak 
occur ing  at  X  . 

Thus  by  varying  m,n  so  that  X  lands  at  or  near  the  w 

c  j 

we  can  detect  all  the  frequencies  provided  the  process  is  purely 


harmonic . 

Again,  as  in  the  previous  case,  in  the  presense  of  noise  this 
procedure  breaks  down  and  the  expected  normalized  HOC  converge  to 


X  . 
c 

2 

Theorem  3.4.  Assume  a  >0.  Then  for  m,n  such  that  X  ■ 

-  c  c 

cos  1 ( (n-m) / (n+m) ) ,  we  have 


11m 

k-*a> 


nZD 

Hk 
N- 1 


V 


Proof .  Define  a  sequence  of  probability  density  functions  on 
[0,n]  . 

..  /,  V  -  [(1-cos  X)®(l+cos  X)n]k  _  , 

(pk{X)  n  - k  *  1 . 

!  [(l-cos  X)®(i+cos  X)n] kdx 
J  0 

Then  for  every  e  >  0  it  is  not  difficult  to  see  that 

*k(X)-+°,  k — »*>, 

J|X-X  | 

'  c 


from  which  follows  that 


(3.3) 


*>k(X  )  cos  ( X  )dX— » 

0 


cos(X  ) 
c 


n-m 

n+m 


Define  now 

h( X )  «  (l-cos  X )m( l+cos  X)n 

and  without  loss  of  generality  assume  0  <  \^  <  n.  Consider  two 

cases.  First,  suppose  w  *  X  Then  obviously  h(w.)  <  h(X  ),  J 

r  c  j  c 

r,  ar.d  (h(w.)/h(X  ))k — >0,  k — ►  <».  Also,  for  sufficiently 
J  c 

small  b  >  0 
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hk(X)dX 


H  »  « 

}  hk(U<«  i,»_,  |>6  *>k(X,dX  ♦  J(l_( 


0  s 


»k(»c) 


Vs6 


•>k<V 


max(hk(X  -6),hk(X  +6)*n 

*  - S- - E -  +  26— >26,  k— »a> 

*» 

and  let  6 — »0.  We  have  shown  in  fact  that 


<pk(wr)— +<K'  k—+®' 


r*ED 


Therefore,  from  (3.2)  cos 


Hk+1^ 


N-l 


coa  x  ,  k — »oo. 
c 


When  * 


Xc>  for  all  j,  then  by  a  similar  argument 


♦  °,  k — » ® 


so  that  from  (3.2)  and  (3.3)  again  coa 
in  general 


nr  ED. 


k+1 


N-l 


cos ( X  ) 
c 


Thus 


uED 


Hk 

H^T 


From  the  last  two  theorems  we  have 

2 

Corollary  3.4.  a ^  *  0  if  and  only  if  for  j 

«ED, 


1, 


H. 


N-l 


Xc,  k— CD, 


3.3.  The  alpha  filter  and  the  case  of  a  single  frequency. 

Although  the  results  of  the  previous  two  subsections  are 

2 

somewhat  pessimistic  when  o  c  >0,  the  fast  conclusion  that  the 

addition  of  noise  makes  the  detection  of  w,  ,  .  .  . ,«  from  HOC 

1  P 
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Impossible  is  far  from  being  true.  Successful  detection  can  be 
achieved  even  in  the  presence  of  appreciable  noise,  provided  care¬ 
fully  designed  filters  are  used  in  generating  useful  HOC  sequences 
A  clue  to  this  effect  is  furnished  by  Theorem  3.4.  The  theorem 
shows  that  in  the  presence  of  noise,  the  normalized  expected  HOC 

converge  to  A  .  The  precise  reason  for  this  fact  is  that  as 
c 

k — ►  ® ,  more  and  more  spectral  weight  is  given  to  \  ,  rendering 

it  dominant.  As  a  result  the  sequence  nED„  /(N-l),  k  -  1,2,..., 

Hk 

is  attracted  to  Ac  and  convergence  occurs.  When  \c  coincides 
with  an  the  resulting  sequence  of  expected  normalized  HOC 

will  converge  to  it.  This  shows  that  by  controlling  and  shifting 
the  spectral  mass  we  can  force  the  sequence  of  normalized  HOC  to 
converge  to  desired  frequencies. 

More  generally,  from  (2.1)  we  obtain  the  basic  spectral  repre 
sentation  for  the  expected  number  of  zero-crossings  (since  Wj  >  0 


(See  Kedem  (1986)  for  a  discussion  and  additional  references  con¬ 
cerning  this  representation. )  From  this  representation  we  see 
that  nED1/(N-l)  is  a  weighted  average  of  the  spectral  support. 
Therefore  nED1/(H-l)  will  change  its  location  with  shifts  in  the 
spectral  weight  dF(<y)  ,  and  so  this  quantity  can  be  "directed"  to 
admit  values  near  or  at  discrete  points  in  the  spectral  support. 

The  point  of  this  discussion  will  now  be  demonstrated  in  the 
special  case  when  p  -  1.  We  will  show  that  <•>  can  be  detected 


by  a  certain  sequence  of  expected  normalized  HOC  that  converge  to 


regardless  of  the  magnitude  of  the  signal  to  noise  ratio. 

Let  p  *  1 , 

■  A  cos(«1t)  +  B  sinfOjt)  ♦  c 
2 

where  A,B  are  Independent  NtO.cr^)  random  variables  and  lnde- 

2 

pendent  of  the  white  noise  {e  ~  tHO.o^).  Define  the 

a-fllter  by 

Y  »  ( 1-a ) Zt  +  aY  ,  t  -  0,±1 . 

and  -1  <  a  s  1.  The  squared  gain  of  this  filter  Is  given  by 

2 

|H  (o)|2  ■=  - - x,  -n  <  «  s  n. 

l-2a  cos  u>  +  a 


and 


n 


cos(o) |Ha(w)  |  dw 


an 


1-a 
1+a ' 


Therefore,  from  (3.4) 

n 


cos 


nED 


H 


N-l 


[  cos(o)|H  («)|2dF(<o) 
J  n  ^ 


r  2 

J  |Ha(«)rdF(») 


n  n 

[  cos(u)|H  (u)|2dFv(u)  +  [  cos (w ) | H  (u)|2dF(v) 
Jq  a  *  Jq  a  c 


n  n 

|  |Ha(«) I  dFx(u)  +  J  |Ha(w) |2dF£ («) 


(3.5 


2  2 

(l-a)aj  a£ 

- *  cos«V  +  n 


l-2a  cos  c^+a 


( 1-a )a\ 


1-2 a  cos  Wj+a' 


1+a 


and  we  see  that  cos 


■nED 


Ha 


N-l 


Is  a  weighted  average  of  cosfu^) 


and  a.  Suppose  s  cos(Wj) .  Then 


-5  cos 


nDH 

r  a  ' 
0 


Tf=T 


S  COS  (<*>J  ) 


and 


nED 


cos 


H 

ao 


nED 


N-l 


5  COS 


H 


N-l 


4  COSfUj ) 
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!.  >>>v>5r>.f  L  _ asa^fe 


or  more  generally 


CIq  4  Oj  S  •  .  .  <  Clj  S  ...  5  C09(Wj). 

Thus  (a^)  a  monotone  Increasing  and  bounded  sequence  which 

converges  to  a  say.  But  then  from  (3.5)  we  have 

2 

1  c  * 

-y  (a  -  cos(«j 


(3.6) 


a  -  cos (Uj ) 


1  •4-Ct 


*  2 
(1-ct  )af 


.2  ' 


*  *  2 

1  -  2a  cos (cj^  )  +  a 


1+a 


Suppose  a  -  cos(v>^)  e  0.  then  dividing  both  sides  of  (3.6)  by 


a  -  coslUj)  leads  to 


1-fa 
1  2 


<  1 


(1-a  )o 


a 


2a  cos  o.  +  a 


*"2  + 


1-fa 


and  hence  to  a  contradiction.  Therefore 


a 


cos (Uj ) 


or 


nED 


V 


N-l  =  W1 


When  aQ  >  cos(w^),  the  sequence  (a^)  is  monotone  decreasing 
and  bounded  and  by  the  same  argument  converges  to  cos(u^ )  from 
above . 

From  the  proof  of  this  theorem  we  can  see  that  any  filter 
with  transfer  function  H  (<•» )  which  depends  on  a  parameter  Q  r- 

\7 

(-1,1)  can  be  used  in  the  detection  of  '•> ^ .  All  that  is  needed 


i s  that 


When  p  =  2  we  can  prove  in  the  sane  way  that  for  large  j , 

nED  /  { N—  1 )  admits  values  between  u.  and  c*>„  . 

H  1  Z 

aj 

Theorem  3.6.  Let  p  *  2  and  choose  an  aQ  e  (-1,1).  Then  regard- 

2  2  1/2 

less  of  the  signal  to  noise  ratio  (cr^-*-a2)  c  ' 


or  equivalently 


— >ct  e  (cos  «2,cos  Uj] 


nED 


H 


N-l 


w  e  [w1  ,a>2  ]  . 


Proof .  As  before  a  is  monotone  and  bounded  and  thus  converges 


J 


to  a ,  where 


cos  w2  <  a 


2  2  2  2 
a1(l-2a  cos  «2+a  )cos  o1+o2(l-2a  cos  )cos  c>2 

5  “  5~  5  5 

ai(l-2a  cos  «2+a  )+a2(l-2a  cos  a>1+a  ) 


cos  w. 


□ 


3.4.  Detection  of  periodicities  by  HOC . 

The  generalization  of  Theorem  3.5  to  the  case  p  2  2  requires 
more  sophisticated  filtering  which  we  shall  not  pursue  here.  In 
practice,  the  generalization  takes  a  somewhat  different  route. 

Since  the  normalized  HOC  tend  to  admit  values  near  or  at  dominant 


frequencies,  the  central  idea  is  to  evaluate  the  periodogram  at 
the  normalized  HOC.  The  combination  of  HOC  and  the  periodogram  in 
this  manner  has  been  reviewed  and  discussed  in  Kedem  (1986). 


4.  A  complete  solution 

Theorem  3.5  shows  that  In  one  special  case  It  Is  possible  to 
determine  a  single  frequency  from  HOC  regardless  of  the  signal  to 
noise  ratio.  In  this  section  we  give  a  general  solution  to  the 
problem  of  determining  p,  , . . . In  the  presence  of  any 
colored  Gaussian  noise  {?  }  with  continuous  density  f(o>).  But 
rather  than  using  HOC  directly  we  use  functions  of  HOC.  Generally 
speaking,  it  Is  sometimes  more  beneficial  to  use  functions  of  HOC 
and  in  particular  the  correlations  {p(k)}  which  by  (2.3)  are 
functions  of  expected  HOC.  Thus,  we  will  show  that  the  oscillation 
in  (Z^)  as  depicted  by  the  expected  HOC,  obtained  by  repeated 
differencing,  determines  the  discrete  frequencies  and  their  number. 
Recall  (2.3)  and  define 

hN(\)  ■  h^XjDj . Cj|_1 ) 

N-l 

■  N"3/4[1  +  2^T  . »n;n)coB(nX)] 

n*=l 

N-l 

*  N_3/4[l+2^  (l-^)Pn(1t> . n*> ;  o )  cos  ( nX  )  ] . 

n-=l 

We  shall  investigate  the  asymptotic  behavior  of  (h  (X))  and  show 
that  this  sequence  of  functions  of  expected  HOC  determines  p  and 

Wj . regardless  of  the  signal  to  noise  ratio  where  the  noise 

1 8  any  colored  noise.  This  is  done  by  showing  that  as  N — ►  <», 
the  sum  represented  by  h  (X)  vanishes  for  X  *  <»>  ,  j  *  l,...,p, 

N  J 

but  diverges  otherwise.  In  this  regard  the  key  idea  is  the  defi¬ 
nition  of  an  Ag-interval. 

Definition  4.1.  Let  (g..(X))  be  a  sequence  of  continuous  func- 

-  ti 


tions  on  [(>,«].  Let  AQ  be  a  positive  number  and  1^  ( o.  )  a 


subinterval  of  [0,n].  We  say  that  I„  -  (a, (i)  is  an  Art-interval 

N  0 


°f  g„(A)  if  the  following  definitions  are  satisfied. 


(a)  gN(X)  <  Aq/  X  e  iN 


(b)  9N<a)  ■  9n(^) 


(c)  The  Lebesgue  measure  of  I„  >  N_1,/2. 

N 


We  will  show  that  asymptotically,  the  number  of  Ag-intervals 


of  h^fA)  is  equal  to  p-1 . 


We  start  off  with  a  lemma  due  to  Wang  (1983).  The  upper 


bound  given  here  is  an  improvement  over  the  one  given  there. 


Lemma  4.1.  Let  (*((n),  n  =  0,±1,...}  be  a  real  valued  stationary 


Gaussian  process  with  mean  zero,  and  square  lntegrable  spectral 


density  function  f(A).  Then  for  any  a  *  o  and  k  such  that 


k  =  2a  +  3/2  we  have  the  inequality 


«  - - 

e|suP  *  ( 4»  +  /S5F)/ 2j”  f2(X)dX. 


Proof .  For  every  fixed  X 


(  n)  naein*  |  2  =  (m)m2a  +  (m)?  (m+1  )ma(m+l  )acos  X  + 


.  .  +  2K  (N)t;  (D^cosCN-lJx 


H-l  N-n 


I  x  K  (m)< (m+n)ma(n+m)acos(nX ) 


n=0  m= 1 


N-l  N-n 


I'l  K  (m)K  ( m+n ) ma ( n+ra ) a | 


n=0  m*l 


N-l  N-n 

5  2^T  !  (mK  (m+n) -R(n)]ma(n+m)a|  + 

n=0  m=l 
N-l  N-n 

(4.1)  +  2^  |  ^R(n)ma(m+n)a( 

n=0  m=l 

where  R(n)  s  (m)?  (m+n)  .  Note  that  the  last  two  expressslons  In 
the  sequence  of  inequalities  are  independent  of  \ .  The  first  of 
these  terms  can  be  simplified  by  introducing  the  lag-process 

Yn(m)  »  <(raK(n+ra)  -  R(n),  m  =  0,±1,..., 

where  EYn(m)  =  0  and  by  the  Gaussian  assumption 

EYn(m) Yn(m+k)  =  ZK  (m)< (m+n)? (m+k)£ (m+k+n)  -  R2(n) 

=  R2 ( n)  +  R2  ( k )  +  R( n-k) R( n+k)  -  R2(n) 

=  R2 (k)  =  R(n-k)R(n+k) 


so  that  (Yn(m))  for  each  fixed  lag  n  is  wide-sense  stationary 
with  mean  zero.  Also,  if  we  define 


Ry  (k)  =  EYn(m)Yn(mtk) ,  k  =  0,±1 . 

n 

then  by  Cauchy-Schwarz  inequality  and  orthogonality  of  sines  and 
cosines 


CD 


£  |RV  (k) (  s  y  ( R2 ( k ) +R2 ( k ) ) 


k=-<x> 


n 


k=-<3D 

oo  jn 

v  2 

a  2  y  R  (k)  < 

k=-oo  "  -n 


f  (1  )d* 


Therefore  (Y  (m))  has  a  continuous  spectral  density  f  (\), 


say,  such  that 


(4.2) 


sup  f  (X)  s  4n  f  (Q)dct), 
X  0  -rr 


uniformly  in  n.  Going  back  to  the  first  of  the  last  two  expres¬ 
sions  in  (4.1)  and  by  invoking  in  succession  first  Cauchy-Schwarz 
inequality  and  then  (4.2)  we  have 


El  y  [g  (m)t  (m+n)  -R{  n)  Jma  ( m+n  )a 


=  E|  ^Yn(m)ma(m+n)a|  s  |e  |  ^  Yn(m)ma(m+n)a  I  2| 


n  N-n 


.....  < 

f  .ST1  a.  .  .a  imx.2,  ..  ...'I 
=  |  I  (m+n)  e  |  fn(X)dX| 


-n  m=l 


w  N-n 


<  |4n  f 2 (X )dX  | y ma(m+n)aelmX I 2dxj 


-w  m*l 


i 

s  {«*j 


f  (X  )dx  •  27i  ■ 


N-K^-N20}1 


/  r  2 

=  2rr  /  2J  f  (X)dX-N 


2a+l/2 


For  the  last  expression  in  (4.1)  we  have 


N-l  N-n 


X'l  R(n)raa(a+n)a|  <  |R(n)  |  ^ma(m+n)a 


n=0  m=l 


Returning  to  (4.1)  equipped  with  these  new  expressions  which  are 
independent  of  \ ,  we  see  that 


Next  we  show  that  hN<\)  tends  to  be  inflated  as  N — ►  ®  in  small 
neighborhoods  of  the  and  vanishes  elsewhere. 

Define  the  set  A  by 

P 

A  =  ft  (*■  e  f  0  ,rr  ]  :  N'3/4  s  |\-<y.|  <  2n-N~3/4}. 
j-1  3 


V 


I 


sup  h  (X)  s  =- 
Xej4  N  RZ 


Toy  52-aj  +  (4n+/4n)/  2j  f2(X)dX  N_ 

L  j=l 


holds . 
Proof . 


Rz(0)hN(\) 


=  m-3/4 


N-l 

1+2  ^  •  *  •  /®n;n)cos(nx  )  r, 

n=l  J 


=  m-3/4 


A 

Rz<0)+2^T  ( l-J)R2(n)cos(nX ) 


=  «“7/4 


^Rz(0)+  2y  ( M-n ) R7 ( n ) cos ( nX  ) 


N  7/4e|  ^ ZnelnX  y  Zne~lnX 


N  7/4E| y Zneln^ I 2. 


Define 


XU)  -  l£zn. 

m=l 


inX  .  2 


sin2 

[*] 

sin2 

ix 

Then 

N 


EiIv 

n=l 


-inX  .  2 


n  N 


rr  N 


'  inw  -inX . 2 
e  e  I  dP 


-n  n*l 


x'"'+|  'Z 


in(<.>-X  )  .  2^  .  , 
e  7  |  f(a>)dw 


—77  n«l 


p  r 

1  V  2* 

2  2.aJ 


in(«  -X ) 


^  in(«.,+X)  1  R 


-inX 


i  ri  ri  #•# 


$1 

I 

1  ® .  tl 


,W7 


& 

t*.‘V 


j.j  tsln  2(wj~x)  8ln  2'<JJ+X) 


Now  observe  that  since  r?(X)  has  a  maximum  at  A  =  0, 


inf  77(A)  -  n(w/M)  »  - -  > 

Ae ( 0 , n/n )  sin  n 


Therefore,  over  the  X -interval  |  A.  — ^  {  <  n/ N 


,-7/4 


V‘>  ‘  R^oT  E 


Iv 


-inX 


N_7/4  1  2 


>  rTToT  2ajr?(x‘wj) 


.  N-7/4  1 _2  4N2 
RJTO)  2  J  ff2 


or,  for  a  fixed  J 


2o2N1/4 


inf  h  (A)  >  -jJ- - 

ia-^k}  "  rz(0) 


To  prove  the  second  part  of  the  theorem  observe  that  for  A  e  d 


and  N  that  satisfies  l.,2.,3.  we  have 


l„-3/4  1,.  .  _  1-3/4 

2N  5  j lXtWj I  5  n  ~  JN 


By  using  Lemma  4.1  with  a  *  0  we  finally  obtain  for  A  e  A 


,-7/4 


hN(X }  =  Rz(0)  E| 


X‘n' 


-inA 


1 


1 


n-7/4 
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and  note  that  the  last  expression  is  independent  of  X  e  4.  r? 

With  the  help  of  this  theorem  we  can  finally  determine  p  and 
the  w's  by  constructing  the  A^-intervals  of  h^(X)  and  letting 
N  increase.  It  Is  easy  to  see  that  for  any  fixed  N,  the  number 
of  A0~intervals  of  h^ ( A )  is  finite  and  that  the  A^-intervals 
are  well  defined.  This  is  the  subject  of  the  next  result. 

First  we  obtain  p. 

Theorem  4.2.  Let  AQ  be  an  arbitrarily  chosen  positive  constant 
and  define 

p„  c  l  +  the  number  of  A  -intervals  of  h„(X). 

N  ON 

Then  regardless  of  the  signal  to  noise  ratio 
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Proof .  By  Theorem  4.1,  when  N  is  sufficiently  large  and  for  1 


such  that 


*  «  [<*>d  +  N  3/4.  w1  +  1-H"3/4], 


1  s  i  s  p-1 


we  have 


hN(\)  <  Aq, 


(wi+i"<t>i>  _  2N  3/4  *  Bln  (wi+rwi)'2N  3/4  >  N  1/2- 

js p-1 

On  the  other  hand  for  large  N 

lim  hN(\)  >  Aq. 

Also,  note  that  the  intervals  where  h^lX)  must  cross  AQ  at 
least  once  are  such  that 


<  (od+N  3/4)  -  (wi+^)  =  N"3/4(1-tiN  1/4)  <  N  1 


-1/2 


0  <  (o1-n/H)  -  (w1-N~3/4)  -  N  3/4(l-nN~1/4)  <  N  1/2. 

It  follows  that  between  and  wi+1'  i  *  l,...,p-l,  there  is 

exactly  one  AQ-interval  of  hN(\),  provided  N  is  large  enough 

Now  consider  the  two  extreme  points  Uj,Up.  Assume  that  0  <  Oj 

•>  <  n .  Then  for  large  N 
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sup  h  (\)  <  A 

-3/4  N  0 
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hN ( X  )  <  A0 • 
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Therefore,  h ..  ' \ )  does  not  crosj  level  AQ  over  [0,w  -N_3/4)  u 

~3/4  * 

(wp+N  <  \  5  n]  and  consequently  hN(\)  has  no  A0-lntervals 

over  [0,u  )  and  over  (o  ,n]  .  Thus  pu-»p,  N — ►  <x> . 

From  the  proof  of  Theorem  4.2.  It  is  clear  that  the  AQ- 

lntervals  of  h^ ( A ) 


1l.tf  S  (bj,N'aJ  +  l,N)  ’  J  *  *'2 . pn-1 


satisfy  for  large  N  and  for  some  a  ,bn  ,  0  s  a,  „  <  w  < 

1 ,N  p.N  1 ,K  1 

bl,N  <  *2 ,  N  <  “2  ‘  b2 ,  M  <  •  •  <  Vl.N  <  ap,N  <  “p  5  bp.N  '  " 
where 


max  {b.  -a .  <  2H 
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We  finally  have 


Theorem  4.3.  Suppose  (a^  ^},{bj  N }  correspond  to  the 
AQ-intervals  of  hN(X).  Let 


Then  for  (3  <  3/4 


"j  '  la3.H'bj,N)- 


11m  H3!  w |  =  0 
N-*x>  J  J 


S?: 


regardless  of  the  signal  to  noise  ratio. 
The  last  three  results  show  that 


1 


Corollary  4.1.  In  the  Gaussian  case,  p,^ . «  are  completely 

determined  by  <EDj}jal  provided  the  spectral  density  of  the 
noise  term  is  square  integrable. 
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